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JLpSTRACT  fConUnue 

■^Models  of  the  electrical  current  collected  by  a planar  aperture  ion  sensor 
generally  assume  that  the  2jr  steradian  field  of  view  is  unobstructed.  This  i^s 
not  always  achieved.  This  report  shows  the  development  of  a model  for  the 
measured  current  based  on  a simple  approximation  for  the  obstructions  that 
are  placed  in  the  field  of  view  of  any  planar  aperture  ion  sensor.  The  details 
of  the  model  given  herein  are  related  to  the  ion  sensor.s  onboard  two  Air  Forc« 
satellites,  where  the  sensors  have  been  recessed  below  the  level  of  the 
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'JO.  Ab-itfact  (Continued) 


■lipai'ei  rafl  .skin.  However,  the  model  can  oa.sily  be  adapted  to  many  other- 
ca.ses  of  ob.st  ruction.s  in  the  field  of  view. 

Kor  Collection  of  data  when  the  ambient  ion  thermal  .speed  is  negligible 
compared  to  the  relative  bulk  flow  (ram)  speed,  the  obstructions  are  (rnly 
important  when  the  flow  vector  from  the  ambient  plasma  to  the  plasma  sensor 
intersects  the  obstructions.  For-  the  collection  of  data  when  the  ion  thermal 
speed  is  comparable  to  the  r-elative  bulk  flow  speed,  'he  obstructions  affect 
the  amount  of  current  collected  by  the  sensors  at  all  orientations  of  the  space- 
craft to  the  bulk  flow-  direction. 
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A Model  for  the  Electrical  Current  Collected 
by  a Planar  Aperture  Ion  Collector  With  a 
Partially  Blocked  Field  of  View 


I.  INTROm!».TIO\ 

Instruments  designed  to  measure  the  bulk  motions  of  environmental  ions  in 
the  upper  atmosphere  were  flown  on  two  Air  Force  scientific  satellites.  Each 
instrument  consisted  of  an  array  of  four  identical  planar  electrostatic  analyzers. 

A planar  electrostatic  analyzer  measures  the  flux  of  positive  ions  through  an 
aperture  from  a field  of  view  of  2 tt  steradians.  The  amount  of  positive  flux  or 
current  measured  is  a function  of  the  analyzer's  efficiency,  its  aperture  area,  the 
ambient  ion  density,  the  mean  thermal  speed  of  the  ions,  and  the  bulk  flow  velocity 
of  the  ions  relative  to  the  plane  of  the  aperture.  The  aperture  planes  of  the  anal- 
yzers or  sensors  in  each  instrument  are  separated  from  one  another  by  approxi- 
mately 40  deg.  If  the  mean  thermal  speed  of  the  ions  and  the  velocity  of  bulk 
plasma  flow  are  significantly  less  than  the  satellite  velocity,  the  difference  in  cur- 
rent flowing  to  the  various  unobstructed  sensor  apertures  can  be  used  to  determine 
the  spacecraft  attitude.  ^ Conversely,  if  the  spacecraft  attitude  is  known,  the 


(Received  for  publication  22  April  1077) 

1.  Sagalyn,  R.C.,  and  Smiddy,  M.  (19G9)  Positive  ion  measurements  of  space- 
craft altitude:  Gemini  X and  XII,  J.  Spacecraft  and  Rockets  6(9):985-991. 
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instrument  can  be  used  to  determine  the  mafjnitude  and  direction  of  bulk  motifins 

i' 

of  the  environmental  ions. 

!i 

The  two  spacecraft,  referred  to  here  as  satellites  A and  H,  were  launched 
into  high  inclination  orbits  with  perigees  of  approximately  250  km.  .Satellite  A 

has  an  apogee  of  1500  km,  and  satellite  B has  an  apogee  of  liOOO  km.  Both  space-  ; 

craft  also  carry  boom-mounted  spherical  electrostatic  analyzers  to  determine  the 
density  and  temperature  of  environmental  electrons.  Also,  onboard  sensors  con- 
tinuously determine  the  spacecraft  attitude. 

Previously  developed  methods  to  analyze  the  data  from  similar  instruments 
have  had  to  be  significantly  modified  for  two  reasons.  First,  the  fields  of  view  of 
the  individual  sensors  are  all  somewhat  less  than  2 ir  steradians.  This  is  due  to  a 
spacecraft  design  compromise  that  led  to  the  placing  of  the  instrument  below  the 

level  of  the  spacecraft  surface.  Second,  over  a considerable  part  of  both  space-  i 

craft  orbits,  the  mean  ion  thermal  speed  cannot  be  taken  as  negligiV)le  compared  I 

to  the  spacecraft  velocity. 

This  report  shows  the  development  of  a procedure  for  analyzing  the  data  from 
the  ion  bulk  motion  instruments.  The  current  to  each  sensor  is  modelled  based  on 
the  obstruction  of  the  field  of  view  and  the  fact  that  the  ion  thermal  .speed  may  be 
comparable  to  the  satellite  velocity.  Although  the  development  is  specific,  it  is 
easily  adapted  to  other,  similar  instruments. 

2.  INSTIU  \11:M  \TI0\ 

The  basic  design  of  the  planar  electrostatic  analyzer  in  the  bulk  motion  instru- 

2 3 4 5 

ments  has  been  given  in  detail  elsewhere.  ’ ' ’ ' Below,  a brief  description  of 
the  instrument  and  its  sensors  is  given  to  aid  the  reader  in  understanding  the  devel- 
opment of  the  modelling  of  the  currents,  which  is  given  in  Section  3. 

The  instrument  with  its  four  sensors  is  shown  in  Figure  1.  Each  sensor 
measures  the  flux  of  positive  ions  through  the  aperture.  All  the  apertures  have 
1 equal  area.  The  normals  of  the  aperture  planes  are  set  approximately  40  deg 


2.  Lai,  S.T.,  .Smiddy,  M. , and  Wildman,  P.J.L.  (1976)  Satellite  sensing  of  low 

energy  plasma  bulk  motion,  Proc.  AIAA  Conf.  on  Aerospace  and  Aeronaut- 
ical Meteo rology  and  Symposium  on  Remote  Sensing  from  Satellites. 
Melbourne.  Florida,  p.  302. 

3.  .Sagalyn,  R.C.,  and  Smiddy,  M.  (1967)  Positive  Ion  .Sensing  .System  for  the 

Measurement  of  Spacecraft  Pitch  and  Yaw.  Air  Force  D-10  Experiment 
Flown  on  Gemini  X and  XII,  technical  memorandum  AFCRL-07-0158. 

4.  Wildman.  P.  J.  L.  (1977)  A low  energy  ion  sensor  for  space  measurements 

with  reduced  photo-sensitivity,  to  appear  in  Space  Sci.  Instrum.  3(No.  4). 

5.  Wildman,  P.J.L.  (1976)  Studies  of  Low  Energy  Plasma  Motion:  Results  and 

a New  Technique,  AFGL-TR-76-0168, 
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upai-t  ffom  one  anotlier.  Two  .sensors  of  eacti  'nstrument  look  into  ttie  spin  plane, 
and  two  sensors  look  out  of  this  plane.  The  satellite's  nominal  spin  axis  is 
designated  as  the  minus  y-axis.  The  spin  plane  nominally  lies  in  the  satellite  x-z 
plane.  The  ty-axis  is  parallel  to  the  satellite's  orbital  momentum  vector 
(Figure  2). 


Figure  2.  Clrientation  of  .Satellite  With  Respect  to  Its  Orbit  Path 
and  the  Earth's  Surface 


On  satellite  A,  there  are  two  instruments,  one  facing  the  +x  direction  of  the 
satellite  and  one  facing  the  -x  direction.  The  four  sensors  in  the  instrument  fac- 
ing the  +x  direction  are  numbered  as  sensors  1 through  4.  Numbers  2 and  4 look 
into  the  spin  plane.  Number  1 looks  40  deg  out  of  the  spin  plane  toward  the  -y 
axis,  and  number  3 looks  40  deg  out  of  the  spin  plane  toward  the  +y  axis.  The 
four  sensors  in  the  instrument  facing  the  -x  direction  are  numbered  as  sensors 
5 through  8.  The  orientation  of  these  sensors  within  the  instrument  is  identical 
to  the  orientation  described  above  for  the  +x  facing  instrument.  Sensor  number  5 
corresponds  to  sensor  number  1,  etc.  On  satellite  B there  is  only  one  instrument; 
it  faces  the  -x  direction  of  the  satellite.  The  sensors  are  numbered  1 through  4 
in  the  same  manner  as  those  in  the  instrument  on  satellite  A. 
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Because  (he  sensors  in  an  instrument  form  a concave  recess  (Figure  1)  and 
because  tlie  instrument  itself  is  recessed  into  the  satellite  surface,  the  aperture.-, 
of  the  sensors  are  below  the  plane  of  the  spacecraft  surface.  The  center  of  the 
aperture  of  the  sensors  looking  into  the  spin  plane  is  4K  mm  below'  the  surface, 
and  the  center  of  those  looking  out  of  the  spin  plane  is  37  mm  below  the  spacecraft 
surface.  Thus,  unlike  planar  sensors  that  are  flush-mounted  with  a spacecraft 
surface,  the  sensors  in  the  ion  bulk  motion  instrument  have  a significant  portion 
of  their  field  obstructed  by  the  spacecraft  surface,  f igures  3 and  4 show  the  fields 
of  view  for  out-of-plane  and  in-plane  sensors,  respectively.  For  cases  where  the 
Mach  number  of  the  satellite's  motion  relative  to  the  ambient  plasma  is  large 
(M  » 1),  the  current  to  a sensor  is  severely  limited  when  the  plasma  flow  vector 
intercepts  a portion  of  the  spacecraft  surface  before  reaching  the  sensor  aperture. 
When  the  flow  vector  does  not  intercept  the  surface,  the  current  is  similar  to  that 
detected  by  an  "unobstructed"  sensor.  For  the  case  when  the  Mach  number  of  the 
satellite  motion  relative  to  the  ambient  plasma  is  near  one  (M  > 1),  the  current 
measured  by  a sensor  is  less  than  would  be  measured  by  an  "unobstructed"  sensor 
for  all  directions  of  the  flow  vector.  The  degree  of  the  reduction  depends  upon  the 
orientation  of  the  obstruction  to  the  flow  direction. 


Figure  3.  Field  of  View  of  Planar  .Sensors  Look- 
ing Out  of  the  .Spin  Plane 
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SENSOR  HORIZON 


Figure  4.  Field  of  View  of  Planar  Sensors  Look- 
ing In  the  Spin  Plane 

1.  C\I,CI  LMION  OK  t;i  UIU.M  TO  \ I’U^TI\L1,^  IlKOI  kKO  I’l.W  \H  SIASOH 

To  model  exactly  the  shape  of  the  obstruction  of  flow  to  the  sensor  (Figures  3 
and  4)  would  require  an  algorithm  that  would  not  be  efficient  for  the  analysis  of 
large  blocks  of  data.  Instead,  the  obstructions  are  approximated  by  a series  of 
wedges  centered  on  the  sensors'  apertures.  The  apertures  are  approximated  by  a 
point  located  at  the  aperture  center.  Figures  5 and  fi  show  the  arrangement  of  the 
wedges,  and  Figures  3 and  4 show  the  image  of  the  wedges  (dotted  lines)  in  the 

fields  of  view.  , 

The  calculations  are  done  in  a sensor-centered  coordinate  system  labeled  as 

the  X . f rn  system.  The  x axis  is  normal  to  the  plane  of  the  sensor.  The  1 
s s’  s •'  s s 

axis  is  perpendicular  to  both  the  x^  axis  and  the  normal  to  the  surface  into  which 
the  instrument  is  mounted.  The  m^  axis  completes  a right-handed  system.  For 
the  sake  of  mathematical  convenience,  the  wedges  are  aligned  with  the  1^  and  m^ 
axis.  Thus,  for  the  sensors  that  look  out  of  the  spin  plane,  it  is  necessary  to 
make  a rotation  of  this  coordinate  system  to  get  the  best  alignment  of  the  wedges 
with  the  actual  obstructions.  For  the  sensors  looking  below  the  spin  plane,  the 
above  coordinate  system  is  rotated  -04  deg  about  the  x^  axis,  and  for  the  sensors 
looking  above  the  spin  plane,  the  rotation  is  -116  deg  about  the  x^  axis. 
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Figure  5.  Orientation  of  the  Two  Wedges  Used  to  Model  the  Obstruction 
for  the  Out-of-Flane  Sensor's  Field  of  View 


Figure  fi.  Orientation  of  the  Two  Wedges  Used  to  Model  the 
the  Obstruction  for  the  In-Plane  Sensor's  Field  of  View 
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The  expressions  for  the  sensor's  axes,  in  spacecraft  coordinates,  arc: 


cos  cos  V 


sin  P cos  Y 


where,  P,  pitch  angle  of  the  sensor's  look  direction  measured  from  the  'x  axis 
of  the  spacecraft,  and  Y . - yaw  angle  of  the  sensor's  look  direction  measured 
from  the  spacecraft  X-Y  plane.  Since  the  instrument  on  satellite  A faces  the 
+x  direction,  and  since  the  second  instrument  on  satellite  A and  the  instrument  on 
satellite  B face  the  -x  direction. 


1,  = (x  X X ,)/magnitude  (x  X x ,)  = 


2 2 2 

-sin  P cos  Y /(cos  Y .sin  P - sin  Y ,) 
s s ' s s 

2 2 2 ' 
sin  Y /(cos  Y sin  P ♦ sin  Y )‘‘ 

S S 3 


for  sensors  facing  the  +x  direction. 


f,  -(x  X X ,)/magnitude  (x  X X ) 


for  sensors  facing  the  -x  direction.  The  third  axis  is: 

/ (cos^  Y sin^  P + sin^  Y )•' 

/ S S s 

. - I 2 2.21 

m - X XI  = ± I -cos  Y sin  Y cos  P /(cos  Y sin  P,  + sin  Y,)^ 

SSSi  S.S  s'  S s s 


is^  Y , cos  P , sin  P , / (cos^  Y . sin^  P , + sin^  Y ,)^ 


For  the  sensors  looking  out  of  the  spin  plane; 


1'  = cos  0*1  +sino*m 
s s s 


m'  -sin  0*1  + cos  o * m„ 
s s .s 


where  o is  either  -64  deg  or  -116  deg. 
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The  general  forinuUi  for  cui'rent  flowing  to  :i  surface  in  a [jlasina  is: 


(Is)  Nef  (T)  (Ir 


(a) 


where  els  surface  element 

N = number  density  cjf  particles 
e electronic  unit  of  charge 

fCr)  particle  distribution  function  in  velocity  space. 

It  is  assumed  that  the  particular  distribution  is  a Maxwellian  distriVmtion; 

a / 9 r 9 9 9 "1 

nr^.i'y,  fj,)  - (a/^)  “ exp  j^-a(r^  * r^,  < , ;)J 

where  a m/2kT  = (4/tc^) 

c average  thermal  speed  of  the  particles 

and  i and  are  the  rectilinear  components  of  the  velocity  of  a particle  in  a 

coordinate  system  at  rest  in  the  gas.  In  order  to  use  Eq.  (5),  it  is  necessary  to 
transform  the  coordinate  system  to  the  satellite  frame.  The  transformation  is: 

»'x  = *’xs  + q cos  9 

Vy  = I'j  + q cos  ? (fi) 

i-z  = I'm  'I 
where 

cos  (9  = X • (-^)/q 
cos  f = 1 • (-q)/q 
cos  rj  - m • (-q)/q 

velocity  of  the  bulk  flow  with  respect  to  the  satellite. 

The  bulk  flow  of  ions  detected  by  the  instrument  is  composed  of  two  compo- 
nents. The  velocity  of  the  satellite  measured  in  geographic  coordinates  (v^)  causes 
an  apparent  flow  of  the  ambient  plasma  in  the  opposite  direction;  this  is  called  the 
ram  velocity.  The  other  component  is  the  true  bulk  flow  of  the  plasma  with  respect 
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to  the  eartii;  tliat  is  tlie  wind.  TJie  total  apparent  flow  velocity  (q)  is  defined  in  a 
non-spin  coordinate  system  such  that  the  x axis  is  parallel  to  the  satellite  velocity 
vector,  the  y axis  is  perpendicvilar  to  both  the  satellite  velocity  and  radius  vectors, 
and  the  z axis  completes  a right-handed  system.  In  this  system,  the  components 
of  the  wind  are  (DX,  DY,  DZ).  The  total  apparent  flow  velocity  can  be  written: 


^satellite  ^wind 


From  the  onboard  altitude  data,  the  pitch  and  yaw  of  the  satellite  v'elocity  is  known 
at  all  times  and  the  bulk  flow  velocity,  expressed  in  satellite  coordinates,  is: 


cos  P cos  V (V  - DX)  + cos  P sin  Y DY  - sin  P DZ 

V V S V V V 


-q 


sin  Y (V  - DX) 
V s 


cos  Y'  DY’ 

V 


(7) 


sin  P cos  Y’  (V  - DX)  < sin  P sin  Y'  DY’  - cos  P DZ 

■ V V S V V V 


It  is  only  necessary  to  substitute  Eqs.  (fi)  and  (7)  into  Eq.  (5)  and  integrate  over 
the  proper  limits. 

In  order  to  maximize  the  speed  with  which  the  integration  can  be  done  with  a 
digital  computer,  the  observed  current  (I  ) is  not  computed  by  direct  integration. 
Instead,  the  current  blocked  by  each  wedge  (1^  and  1^)  is  subtracted  from  the  total 
current  flowing  to  an  unobstructed  sensor  (1^).  Note  that  in  the  case  of  the  sen- 
sors looking  out  of  the  spin  plane,  the  current  from  the  region  where  the  two 
wedges  overlap  is  subtracted  twice.  This  error  is  considered  to  be  negligible, 
since  the  area  of  overlap  is  small  and  is  always  separated  from  the  bulk  flow 
‘direction  by  a large  angle.  Thus,  the  calculated  current  is: 


Is  = Iq  - h - l2 


(8) 


where 


4 

k 


Iq  = 1/2  NeAq  cos  fl[l  + erf(w)  + exp  (-w  )/w\^7r] 
A area  of  the  aperture 
w = Qv'  a cos  0 


1 

r. 


ir> 


1 


For  the  out -of-spin-plane  sensors. 


\ e A (a  ' rr) 


-Of  -i^cotOj 


di  f d,.  (-1  ) 
ni  J X X 


2 2 
exp  (-a{(i  j ■-  q oos  f)  ‘ ^ q cos  n)~  ^ (r^  ‘ q cos  n)  ') 

which  integrates  into 

- 1/2  \eA  (a/")-  {\  qcoshfl  -*  erf(w)  * erp  ( -w^) / ws^tJ 


Of 

-(a/T)^  / I erf  a (v,  cot  o,  ± cos  rj)]  exp  (-a  (v  - q cos  h)“|  dv  ) . 

^ Z Z i Z 7. 

o 

Due  to  the  placement  of  the  wedges,  the  upper  sign  is  used  for  the  sensors  looking 
below  the  spin  plane,  and  the  lower  sign  is  used  for  the  sensor  looking  above  the 
spin  plane.  The  current  is  calculated  in  a similar  manner  to  the  calculation  of 
Ij.  The  result  is  that  the  total  calculated  current  into  a sensor  looking  out  of  the 
spin  plane  is; 

Of 

= 1/2  NeA  {ajir)^  J" {erfK'a  (r^  cot  (<)j  ± q cos  n)| 
o 

r 2 

* erf  [v  a (v ^ cot  <^2  ' 1 ' ^'''7.  ' 

For  the  in-plane  sensors,  the  above  process  is  repeated  with  different  integration 
limits.  The  result  is: 

00 

= 1/2  NeA  (a/TT)^  J"  ferfjv'a  (t/^  cot  + q cos  n)| 
o 

2 

+ erf  (v/a  (v ^ cot  <ii^  - q cos  f )|  } exp  [-a  (r^  - q cos  f))  | dr,,  . (10) 

The  integration  in  Eqs.  (9)  and  (10)  cannot  be  done  analytically,  but  can  be  calcu- 
lated with  good  accuracy  with  a 10-point  Gaussian  quadrature. 
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\\  itli  the  above  equations,  there  ar  e five  unknowns  tc.  tie  iletermineti;  tlie  wind 
velocity  components  (1)X,  O'!  , and  DZ),  tlie  density  (\),  and  t|ie  mean  thermal 
speed  (c  (4  ^a)-).  I'r'om  any  one  data  prdnt,  only  four  simultaneous  currents  are 
available  from  the  sensor  array  near  the  ram  directions.  T hu->,  at  best,  a solution 
is  obtainable  for  only  four'  of  the  five  unknowns.  ( ine  of  the  plasma  parariieter.s 
must  be  set  to  a value  obtaineii  from  either  other  data  or  an  a pidori  assumption. 

If  a series  of  data  points  is  used,  it  is  pos->ible,  in  jirinciple,  to  determine 
all  five  unknown  plasma  parameters.  However,  it  is  not  always  possible  to  deter- 
mine all  five  unknown  parameters  with  an  acceptable  degree  of  uncertainty.  1 here 
are  reasons  for  this  limitation  to  the  data  analysis.  First,  the  data  contains  noise 
from  several  sources.  The  use  of  several  consecutive  data  points  contributes 
more  to  impi-oving  the  signal  to  noise  ratio  than  it  does  to  contributing  information 
not  available  from  a single  data  point.  Secondly,  no  meth<»d  has  been  found  for 
determining  the  plasma  parameters  directly  frcmi  the  data.  The  set  of  equations 
represented  by  (9)  and  (10)  yields  the  expected  data  as  a function  of  the  plasma 
parameters  instead  of  vice  versa.  The  parameters  are  found  b%  a search  of  the 
manifold  represented  by  them  for  the  ''solution”  that  yields  a minimum  in  the  dif- 
ference between  the  expected  and  the  real  currents.  While  there  may  be  a "best" 
solution,  there  may  also  be  a number  of  "good”  solutions  for  widely  separated 
values  of  the  parameters.  In  short,  the  best  available  algorithm  for  finding  the 
plasma  parameters  can  yield  false  answers,  especially  if  the  relationships  between 
the  data  and  the  parameters  are  not  clearly  stated. 

In  order  to  choose  which  four  of  the  five  plasma  parameters  should  be  deter- 
mined from  the  data,  it  is  necessary  to  consider  the  relationship  between  the  model 
current  and  the  various  parameters.  Changes  in  "DZ”  tend  to  shiF  the  model  cur- 
rent curves  on  the  time  axis  with  minimal  change  in  the  height  or  the  shape  of  the 
curves,  ('hanges  in  ' DV”  tend  to  change  the  magnitude  of  the  current  in  the  out- 
of-plane  censor,  without  significantly  changing  the  shape  of  any  of  the  curves  and 
with  a small  change  in  the  magnitude  of  the  current  to  the  in-plane  sensors. 

Changes  in  ' N”  tend  to  change  the  magnitude  of  all  currents  equally,  without  sig- 
nificantly changing  the  slope  of  the  current  curves.  Changes  in  "C”  tend  to  affect 
the  shape  of  the  current  curves  and,  to  a lesser  degree,  the  magnitude  of  the 
current,  without  significantly  changing  the  location  of  the  curves  on  the  time  axis. 
Changes  in  ”DX'  'end  to  affect  the  shape  of  the  current  versus  time  curves  and 
the  magnitude  of  the  current.  Thus,  changes  in  "N"'  and  "C”  can  be  interpreted  as 
changes  in  "DX”  and  vice  versa. 
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In  pracuco,  if  nu  additional  a priori  inrorniati(;n  is  available,  it  is  best  to  fix 
the  value  id'  ' DN"  (to  a value  id'  zero)  and  then  solve  for  the  other  parameters.  If 
simultaneous  data  from  the  spherical  electrostatic  analyzers  onboard  both  vehicles 
can  ^>ive  the  electron  density  .vith  an  accuracy  of  1 0 percent,  the  parameter 
can  be  set  equal  to  the  electron  density,  allowing  ”DX"  to  be  solved  for  from  the 
ion  data,  in  the  case  where  a retarding  potential  analyzer  (HPA)  is  available,  the 
"DX"  IS  determined  from  the  HPA  data. 

■).  ( list  1.1  'KIN' 

In  this  report,  a formula  has  been  developed  to  model  the  current  to  a planar 
electrostatic  ion  sensor  with  a partially  blocked  field  of  view.  The  use  of  the 
fortnula  is  shown  in  Figure  7.  The  data  were  obtained  from  satellite  A when  the 
Mach  number  of  the  satellite  relative  to  the  ambient  plasma  was  approximately 
four.  Kach  second  of  data  repreoents  approximately  19  deg  rotation  of  the  satellite 
about  its  spin  axis,  with  the  ram  direction  passing  through  the  zenith  of  the  sen- 
sor's field  of  view.  Also  shown  is  the  current  curve  given  by  Eq.  (10)  and  the  cur- 
rent curve  given  bv  L,  of  Eq.  (8)  when  qv  a » 1,  that  is,  the  "cosine”  formula  used 
2 ^ 

by  I.ai  et  al.  It  is  obvious  that  the  cosine  formula  is  a poor  model  for  the  data. 

The  complete  unobstructed  sensor  formula  represented  by  1^  of  Eq.  (8)  will  give 
almost  identical  representations  for  the  data  on  the  left-side  and  center  of  Figure  7, 
but  will  give  a value  much  larger  than  the  data  on  the  right-side  where  the  current 
decreases  rapidly  due  to  the  obstruction  of  the  field  of  view.  As  the  Mach  number 
approaches  one,  the  difference  between  the  current  calculated  by  the  formula  for 
the  unobstructed  sensor  and  the  formula  for  the  obstructed  sensor  will  become 
significant  at  all  parts  of  the  current  versus  time  curves. 
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Figure  7.  Comparison  of  the  Best  Fit  Possible  With  the  Unobstructed  Flow 
Formula  and  the  Formula  Obtained  Herein 
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